Identification of a Scalar Glueball 
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We have performed a coupled channel study of the meson-meson S-waves involving isospins (J) 
0, 1/2 and 3/2 up to 2 GeV. For the first time the channels tttv, KK, 7777, aa, gr/' , g'g' , pp, uuj, 
uj(j>, (jxf>, a\TT and n*n are considered. All the resonances with masses below 2 GeV for I — and 
1/2 are generated by the approach. We identify the /o(1710) and a pole at 1.6 GeV, which is an 
important contribution to the /o(1500), as glueballs. This is based on an accurate agreement of 
our results with predictions of lattice QCD and the chiral suppression of the coupling of a scalar 
glueball to qq. Another nearby pole, mainly corresponding to the /o(1370), is a pure octet state 
not mixed with the glueball. 



PACS numbers: 11.80.Gw, 12.39.Fe, 12.39.Mk 

1. QCD, the present theory of strong interactions, 
is a non-abelian Yang Mills theory so that gluons carry 
colour charge and interact between them. It is gener- 
ally believed that QCD predicts the existence of mesons 
without valence quarks, the so called glueballs. Its con- 
firmation in the spectrum of strong interactions is then 
at the heart of the theory. In quenched lattice QCD the 
lightest glueball has the quantum numbers of the vac- 
uum, J PC = 0++, with a mass of (1.66 ±0.05) GeV [1]. 
Experimentally the closest ++ scalar resonances to this 
energy range are the /o(1500) and /o(1710) [2]. Some 
references favour the /o(1500) as the lightest scalar glue- 
ball [3], while others do so for the /o(1710) [4,5]. 

We analyze the 1 = meson-meson S-wave in terms 
of 13 coupled channels, ttit(1), KK{2), 7777 ( 3 ) , aa{4), 
777/(5), 7/7/(6), pp(7), lulo(8), K*K*(9), c^(10) 00(11), 
ai(1260)7r(12) and tt* (1300)tt(13). The number la- 
belling each state is given between brackets to the right. 
The multipion states, which play an increasing role for 
energies above ~ 1.2 GeV, are mimicked through the 
aa, pp and ujlu channels. It is worth stressing that 
our approach is the first one with such a large num- 
ber of channels and that a similar scheme could also 
be applied to other controversial meson-meson partial 
waves. In addition, we study simultaneously the S-wave 
of R-TT+ -> R-TT+ (involving 7=1/2 and 3/2) with 
the coupled channels Kir, Kg and Kg' . 

2. Let T^j be the i j S-wave amplitude with 
isospin / and i,j = l...n, with n the number of chan- 
nels. We use the master formula = [I + N ■ g] 1 • 
N , where N is the symmetric matrix of interaction ker- 
nels and g is a diagonal matrix of elements gi(s). The 
function gi(s) is calculated from kinematics in terms of 
a once subtracted dispersion relation and a subtraction 
constant dj [6]. Since SU(3) breaking is milder in the 
vector sector we take 0,7 = ag = ag = aio = an- The 
rest of subtraction constants are fitted to data. The ma- 
trix elements Nij consist of the sum of two tree level 
contributions. The first is a contact interaction calcu- 



lated from the lowest order Chiral Perturbation Theory 
Lagrangian, €2- The second is due to the exchange of 
bare resonances in the s— channel with the couplings 
calculated from the lowest order chiral Lagrangian in- 
cluding an octet and singlet of ++ resonances, £5 [7]. 
Explicit expressions of Ni j can be found in ref. [6] for 
the simplified case of three channels without includ- 
ing the 771 field. We extend these Lagrangians from 
SU(3) to U(3) as the 771 field is needed to deal with 
the 77 and 77' mesons, similarly as in ref. [8]. The ma- 
trix $ = J2i=i 4>iK/V2 + 'r]i/v / 3 incorporates in a stan- 
dard way the nonet of the lightest pseudoscalars. We 
also employ the matrix U = exp(iv / 2 < i ) //) and the co- 
variant derivative D^U — d^U — ir^U + iUi^, with 
/ the pion decay constant in the chiral limit fixed to 
./tt = 92.4 MeV. The classical left and right external 
fields, and respectively, are necessary to gauge 
the global chiral symmetry to a local one [7] . The field 
= (?> 

since it is identified with A 
of the lightest 1~~ vector resonances and A is a con- 
stant, with A = 4.3 from the width p — > 7T7t. The cou- 
plings of the vector-vector states to the pseudoscalar- 
pseudoscalar and aa ones are then determined by min- 
imal coupling [9]. Our fits require a singlet and two 
octets of bare resonances. The two octets were already 
considered in ref. [8] in the study of K~ir + — > K~ir + . 
We fix the parameters of the first octet, mass and cou- 
pling constants, to those in ref. [8], Mg = 1.29 GeV, 
26 MeV. The bare mass of the second octet 

1.90 GeV. We 



£fj,)/2 plays a special role in our approach 
where is the nonet 



(2) 



is fixed from the same reference, M t 
are then left with three parameters for the singlet, Mi 



Cm , and two for the second octet, c? 



and c~m ■ 



It 



results from our fits that M\ < 0.9 GeV. 

Concerning the aa channel we follow a novel method 
to calculate its transition amplitudes, A^ 4, without in- 
cluding any new free parameter. This can be done be- 
cause the a corresponds to a pole due to the interac- 
tions between two pions in the 1 = S-wave, (7r7r)o 
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[10]. For the interaction kernel TVj 4 one starts by cal- 
culating from the Lagrangians £2 and Cs the tree level 
amplitude T^ s for i — > (7T7r)o(7r7r)o. To take into ac- 
count the pion final state interactions, T^ s is multi- 
plied by the factor JlfeLi ^/D(sk), with m the number 
of cr's in the scattering process (2 or 4) and Sk the to- 
tal centre of mass (CM) energy squared of the k t h pair. 
We use here that the rescattering of two 1 = S-wave 
pions from a production kernel is given by the factor 
l/D = 1/(1 + V 2 gi), with V 2 = (a - m 2 /2)// 2 [10]. To 
isolate TV,-^ one takes the limit (for definiteness i 7^ 4) 
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Where the subscript 77 indicates that the correspond- 
ing function is calculated on the second Ricmann sheet 
(with the sign reversed in the definition of the pion 
three-momentum), s CT is the a pole position and g,,^ 
is its coupling to irir. Performing the Laurent ex- 
pansion around s a of \/Dn(s) = ao/(s — s a ) + ... 
the evaluation of Ni t i from eq.(l) requires the ratio 
(ao/5<T7T7r) 2 . Since gi.n(s a ) = -/ 2 /(s CT - ml/2) at s CT , 
where l+V^ffi,// = 0, and taking Tn ~ V2/ (l+V~2giji), 
appropriate for these energies [10], then {a^/ g^^) 2 = 
/2/(1 _ ^£l| Sg ^-^/ 2 ) 2 ) ~ f\ In this way, 7V M = 

2+S f 2 j ^ 4 ; ^^A AT. . — f 4 



rp-J + O f2 

1 iA J 



and 7V 4 ,4 = TiY p. Using 7V i>4 evalu- 



ated with Sk = So- violates unitarity because s a is com- 
plex and Ni A must be real. Instead, we interpret the 
width of the a resonance as a Lorentzian mass distribu- 
tion around its nominal mass value ~ 450 MeV with a 
width ~ 500 MeV. In this way the a masses (^/sk) used 
to calculate the functions Ni t 4 and 54 are folded with 
the previous mass distribution. Similarly, for the pp 
state gi is also convoluted with a p mass distribution. 

3. We fit our 12 free parameters to 370 data 
points from threshold up to 2 GeV. The data com- 
prise the 1 = S-wave 7777 phase shifts 6®, the elasticity 
t]q = |5i,i|, the 1 = S-wave nn — ► KK phase shifts 
8\ t 2 and modulus |Si,2|, the S-wave contribution to the 
7T7T — > 7777, 7777' event distributions and the phase (0) and 
modulus (A) of the K~ir + — > K~tt + amplitude from 



the LASS data. 

Si.i — fiin ~t~ 2i A 



The S-matrix element Si 



a) 





1000 1500 
(MoV) 




1^1 



1500 
(MeV) 



|0.1 
< 





77-77- —1- 7)7] 




S-wave 




1) A, 




















H .j is given by 

- "ij t "v^^.j v^"' wnere ft = qi/Sn^/s and 
<7i is the CM three-momentum for channel i. In order, 
these data are shown on the first eight panels of Fig.l 
from top to bottom and left to right. For y/s < mx 
in the <5q(s) panel we have the inset showing in detail 
the precise data from if e4 decays. The reproduction of 
the data is fair, as shown in the figure. The dashed 
lines on the first eight panels include the a\ir and tt*tt 
states, while the solid ones do not. The similarity be- 
tween both curves indicates that these channels give 



s/S (MeV) y/s (MeV) 

FIG. 1. Fit to experimental data. More details are given 
in the text. 



small contributions. The width of the band represents 
our systematic uncertainties at the level of two standard 
deviations, n a = Ax 2 /(2% 2 ) 1 ^ 2 [11]. Compared with 
other works [12-14] we determine the interaction ker- 
nels from standard chiral Lagrangians, avoid ad-hoc pa- 
rameterizations, include many more channels and fewer 
free parameters are used. For / = 1/2 the k pole is 
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located at (708 ± 6 - i 313 ± 10) MeV, the X o *(1430) 
at (1435 ± 6 - i 142 ± 8) MeV and the X o *(1950) at 
(1750±20-i 150±20) MeV, similarly to ref. [8]. For / = 
one has the /o(600) or a at (456 ± 6 - i 241 ± 7) MeV 
and the / (980) at (983 ± 4 - i 25 ± 3) MeV. There are 
poles at (1690 ± 20 - i 110 ± 20) MeV, corresponding 
to the / (1710), and at (1810 ± 15 - i 190 ± 20) MeV, 
with mass and width in agreement with those reported 
for the / (1790) by BESII. In the PDG [2] the width 
for the / (1710) is 137 ± 8 MeV, much smaller than 
220 ± 40 MeV from the given pole position. How- 
ever, we have checked that on the real axis the value 
of the width corresponding to the half-maximum for 
the partial waves with prominent /o(1710) peaks is just 
160 MeV [15]. This reduction is due to the opening 
of several channels along the resonance region and the 
agreement with the PDG is restored. The other poles 
at (1466 ± 15 - i 158 ± 12) MeV and (1602 ± 15 - iU ± 
15) MeV, connected with the .f (1370) and / (1500), are 
referred in the following as and f^, respectively. De- 
spite that we have included only three bare resonances 
in / = we have generated six. The poles are located 
on the unphysical Riemann sheets that connect contin- 
uously with the physical one for some interval along the 
real s— axis. Note that the pole /(f does not influence 
the physical axis beyond the 7777' threshold at 1505 MeV, 
since above this energy a different Riemann sheet is the 
one that matches with the physical s— axis. This effect 
typically gives rise to a pronounced signal at the 7777' 
threshold and this is the reason for the /o(1500) mass, 
(1505±6) MeV [2]. If a physical amplitude is dominated 
by the f$ pole, then its peak at 1505 MeV has an effec- 
tive width larger than the one from the pole position, 
88 MeV. This is so because given a Breit-Wigner located 
at the position of the f$ pole the energy interval below 
1.5 GeV at which half the value of the modulus squared 
at 1.5 GeV is reached is 8 = 1.2T = 105 MeV, the width 
of the / (1500) [2]. The / (1370) is mainly given by 
the Jq pole, though its precise shape is sensitive to 
for those channels that couple strongly with the latter. 
In Fig.l we also show in the last two rows data from 
pp inelastic scattering at 450 GeV/c and pp annihila- 
tion by the WA102 and Crystal Barrel (CBC) Collab- 
orations, respectively. We have fitted the WA102 data 
using a coherent sum of Breit-Wigner functions and a 
non-resonant term, similarly as done by the WA102 Col- 
laboration [16]: 

i) sTs <m 7} + m v > , A = {a, /o(980), / L , / *} , 

jeA 1 J J 

ii) sfi > m„ + m„. , B = {a, /„(980), /o(1710), A(l"90)( , 



NR{^s) t = a(VI -m k - m^e"^" 5 * , (2) 

where a,j and 9j are the amplitude and the phase of the 
production vertex of the jth resonance, Mj, Tj and g^i 
are, respectively, the mass, width and the coupling to 
channel i of the same resonance. The latter is deter- 
mined from the residue of the partial waves at the pole 
position. In addition, + mi is the threshold for the 
channel i and a, f3, 7, S are real parameters. The form of 
the non-resonant term is taken from the WA102 Collab- 
oration [16] . The constant ri is fixed so as the amplitude 
A(y/s)i is continuous at w m i = m v + m' v . As explained 
above, once the 7777' threshold is crossed over one has to 
consider other Riemann sheets which do not have the 
/ L and /,f poles but the /o(1710) and /o(1790) ones. 
Above Wr)r)' the a and /o(980) give tiny contributions. 
Tj in eq.(2) is the largest between its value from the pole 
position and the one calculated by summing the partial 
decay widths T Jvl = 6(v r s-m k -mi)\ i \g j ..i\ 2 q i /(8irMj), 
with Xi = 1/2 for identical particles. Eq.(2) incorpo- 
rates important new facts compared to the analyses of 
the WA102 Collaboration. First, the pole positions for 
the different resonances are those already determined 
from our study of the scattering data on the first 8 pan- 
els of Fig.l. Let us stress that these observables only 
involve two particles in the final state and their analysis 
is theoretically cleaner. Second, the couplings gj-i are 
similarly fixed. Third, the a,j and 0j parameters are the 
same for all the WA102 reactions, that are fitted simul- 
taneously. For the Crystal Barrel data on pp annihila- 
tion we also use eq.(2) but without NR(y/s). A good 
reproduction of the data results. In pp — > 7r°?7?7 one ob- 
serves a broad bump for the / (1370) and a prominent 
peak for the /o(1500), that also gives strong signals in 
the WA102 data. Other peaks are observed for the a, 
/o(980) and / (1710). The latter is important for the 
the shoulder in pp — > ppr\r\ above 1.5 GeV. 
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/o(1370) 
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/o(1710) 




3.59 ±0.16 


1.30 ±0.22 


1.21 ±0.16 


\9k«k« 


2.23 ±0.18 


2.06 ±0.17 


2.0 ±0.3 


\dvv\ 


1.7 ±0.3 


3.78 ±0.26 


3.3 ±0.8 


Iswl 


4.0 ±0.3 


4.99 ±0.24 


5.1 ±0.8 


Iswl 


3.7 ±0.4 


8.3 ±0.6 


11.7 ± 1.6 



TABLE I. Couplings of the /o(1370), /(f and /o(1710). 



4. In table I we give the couplings of the f$ (iden- 
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tified as the / (1370)), f£ and / (1710) poles to the 
two pseudoscalar channels. We observe that the cou- 
plings of the Jq and /o(1710) are quite similar. This 
is so because the two poles coalesce in the same one 
when moving continuously from the sheet of one of 
them to the one of the other. They correspond to the 
same underlying resonance, but split in two due to the 
interaction in coupled channels. From the couplings 
of the fo(1710) one can calculate the branching ratios 
T(KK)/T total = 0.36±0.12(0.38t°;°j»), T( m )/T total = 
0.22 ± 0.12(0.18±g;51), and T{tttt)/T{KK) = 0.32 ± 
0.14(< 0.11), where the values of the PDG are given 
between brackets. The values are compatible within 
one sigma. We also obtain that the /o(1790) has a small 
KK coupling, and this is a major difference with respect 
to the / (1710) as stressed by BESII. The couplings of 
the /q'-/o(1370) in table I correspond to the pure 1 = 
octet member (uu + dd — 2ss)/\/6 because they are very 
close to the tree level ones 1^+^- | = 3.9, \gx°R° I = 2-3, 
\g V ri\ — 1-4, \g V r]>\ — 3.7, \g V 'n'\ — 3.8 GeV calculated 
from the Lagrangian Cs [7], with c^m and 
given above. We have also checked that this is the case 
for the -?sTq(1430) resonance which is the 1=1/2 mem- 
ber of the same octet. It follows then that the first 
octet is a pure one without mixing with the nearby /<f 
and /o(1710). The / L -/o(1370) couplings imply a large 
width to tttt with r(/ (1370) -> 47r)/r(/ (1370) -» 
7T7r) = 0.30 ± 0.12, in good agreement with the inter- 
val 0.10-0.25 given in the recent ref. [17]. Let us see 
that the pattern of sizes of the couplings of the and 
/o(1710) corresponds to the chiral suppression of the 
coupling of a scalar glueball, Go, to qq [5]. Accord- 
ing to ref. [5] this coupling is proportional to the quark 
mass, which then implies a strong suppression in the 
production of uu and dd relative to ss from Go- With 
a pseudoscalar mixing angle sin (3 = —1/3 one has that 
7] = -7] s /V3 + Vuy/2j3 and rf = Vs^/V^ + Vu/VS with 
rj s = ss and r\ u = (uu + dd)/V2. Denoting by g ss the 
production of T] s r]s, g sn that of r\ s ri u and g nn for r\ u r\ u , 

gn'v' = 2g ss /3 + g nn /3 + 2V2g„ s /3 , 
g m > = -V2g ss /3 + V~2g nn /3 + g ns /3 , 
g nv = 9ss/3 + 2g nn /3 - 2V2g„ s /3 . (3) 

If the chiral suppression of ref. [5] operates then \g ss \ 3> 
\g nn \- This together with the OZI rule suppress the cou- 
pling g ns . Taking e.g. the couplings of f$ one obtains 
g ss = 11.5 ± 0.5, g ns = —0.2 and g nn = —1.4 GeV, and 
the strong suppression is clear. We now consider the 
KK coupling. A K° in terms of valence quarks corre- 
sponds to y^.-_ 1 siu 1 summing over the colour in- 
dices, and analogously for the K°. The production of a 



colour singlet ss from the K°K° requires then the com- 
bination SiS 1 = 5jss/3 + (sis 3 — Sfss/3), and similarly 
for UjU 1 . As the production occurs from the colour sin- 
glet ss source, only the configuration ss uu contributes, 
picking up a suppression factor of 1/3. In addition, the 
coupling g ss has an extra factor 2 compared to that 
of a ssuu, because the former contains two ss. One 
then expects that the coupling to K°K° has the ab- 
solute value g ss /6. For the and / (1710) it results 
\g K oxo \ — 2 GeV, in good agreement with table I. An- 
other resonance with a known enhanced coupling to ss 
is the /o(980). However, the sizes of its couplings to 7777, 
i]r]' and rj'i]' follow the opposite order to the /o(1790) 
and /(f cases and all of them are much smaller than 
the coupling to KK. Note that quenched lattice QCD 
[4] establishes that the couplings of the lightest scalar 
glueball to pseudoscalar pairs in the SU(3) limit scales 
as the quark mass, in support of the chiral suppression 
mechanism of ref. [5] , that we also observe as discussed 
above. This mechanism also implies that the glueball 
should remain unmixed. This accurately fits with our 
previous result that both the /(f and /o(1710) do not 
mix with the nearby f$. In addition, the masses of the 
fo^ and /o(1710) poles are in excellent agreement with 
the quenched latticed QCD prediction for the mass of 
the lightest glueball, (1.66 ±0.05) GeV. 

5. In summary, we have presented a coupled channel 
study of the I = 0, 1/2 meson- meson S-waves from 
tttt threshold up to 2 GeV with 13 coupled channels. 
All the 1 = and 1/2 ++ resonances with masses 
below 2 GeV have been generated. The /o(1710) and a 
pole at 1.6 GeV, which is an important contribution to 
the /o(1500), are identified as glueballs. Another pole 
at (1.466 — i 0.157) GeV, mainly corresponding to the 
/o(1370), is shown to be a pure octet member. 
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